前言(扯淡)
数论分块是莫比乌斯反演和某些数学题有时候会用到的一个玩意,也是属于那种知道就知道不知道就不知道的内容,现场推如果不知道具体形式感觉还是很难想出来的.这篇博客就尝试用肯定很不严谨的态度尝试"证明"一下相关的内容,实际上就是能在不借助其他信息的前提下能独立推出正确的结论这个目的.
数论分块
考虑求下面这个式子:
∑ i = 1 n ⌊ n i ⌋ \sum\limits_{i=1}^n \left\lfloor\frac{n}{i}\right\rfloor
i = 1 ∑ n ⌊ i n ⌋
数论分块的重要依据就是⌊ n i ⌋ \left\lfloor\frac{n}{i}\right\rfloor ⌊ i n ⌋ 这n n n 个值排成一列,呈单调下降趋势,且是一段一段的.因此可以推出每一段的左右区间各自在哪,并借助类似前缀和的想法快速求和.由于最多只有n \sqrt{n} n 级别的不同的数,因此只要求和部分是O ( 1 ) O(1) O ( 1 ) 的,整个求和部分就是O ( n ) O(\sqrt{n}) O ( n ) 的.
复杂度证明
整个序列的不同的数的个数是n \sqrt{n} n 级别的
对于i i i 而言,有两种取值:
①i ≤ ⌊ n ⌋ i \leq \left\lfloor\sqrt{n} \right\rfloor i ≤ ⌊ n ⌋ 的⌊ n i ⌋ \left\lfloor\frac{n}{i}\right\rfloor ⌊ i n ⌋ 有⌊ n ⌋ \left\lfloor\sqrt{n} \right\rfloor ⌊ n ⌋ 种取值.
②i > ⌊ n ⌋ i > \left\lfloor\sqrt{n} \right\rfloor i > ⌊ n ⌋ 的⌊ n i ⌋ \left\lfloor\frac{n}{i}\right\rfloor ⌊ i n ⌋ 由于⌊ n i ⌋ ≤ ⌊ n ⌋ \left\lfloor\frac{n}{i}\right\rfloor \leq \left\lfloor\sqrt{n} \right\rfloor ⌊ i n ⌋ ≤ ⌊ n ⌋ 因此也只有n \sqrt{n} n 种取值.
综上,分段是n \sqrt{n} n 级的.
假定现在知道区间左端点是l l l ,要求右端点r r r .也就是要找到一个最大的r r r 满足⌊ n l ⌋ = ⌊ n r ⌋ \left\lfloor\frac{n}{l}\right\rfloor = \left\lfloor\frac{n}{r}\right\rfloor ⌊ l n ⌋ = ⌊ r n ⌋ .
r r r 的取值结论是r = ⌊ n ⌊ n l ⌋ ⌋ r = \left\lfloor\frac{n}{\left\lfloor\frac{n}{l}\right\rfloor}\right\rfloor r = ⌊ ⌊ l n ⌋ n ⌋
取值证明
证明他是正确的还是比较容易的.
首先对于这段区间而言,其中每一个值都是⌊ n l ⌋ \left\lfloor\frac{n}{l}\right\rfloor ⌊ l n ⌋
r = ⌊ n ⌊ n l ⌋ ⌋ ≤ n ⌊ n l ⌋ < r + 1 r = \left\lfloor\frac{n}{\left\lfloor\frac{n}{l}\right\rfloor}\right\rfloor \leq \frac{n}{\left\lfloor\frac{n}{l}\right\rfloor} < r + 1 r = ⌊ ⌊ l n ⌋ n ⌋ ≤ ⌊ l n ⌋ n < r + 1
换一下位置可以得到:
{ ⌊ n l ⌋ ≤ n r n r + 1 < ⌊ n l ⌋ \begin{cases}\left\lfloor\frac{n}{l}\right\rfloor \leq \frac{n}{r}\\\frac{n}{r+1} < \left\lfloor\frac{n}{l}\right\rfloor\end{cases} { ⌊ l n ⌋ ≤ r n r + 1 n < ⌊ l n ⌋
这个结论即证明了r r r 就是右端点,右移一位就不满足.
对于初始的l l l 来说,他肯定起始于数列的开头,之后每一次都可以往后推r r r 已经更新l l l 的位置.便可以解决开头的问题了.
模板代码
int f(int n)
{
int res = 0;
for(int l = 1,r;l <= n;l = r + 1)
{
r = n / (n / l);
res += (r - l + 1) * (n / l);
}
return res;
}
例题① [CQOI2007]余数求和
原题链接:luogu 2261
题目大意:
给定n , k n,k n , k ,计算:
G ( n , k ) = ∑ i = 1 n k % i G(n,k) = \sum\limits_{i=1}^n k\%i
G ( n , k ) = i = 1 ∑ n k % i
思路
考虑变换形式:对于k % i k \% i k % i 来说,他等价于k − ⌊ k i ⌋ ∗ i k - \left\lfloor\frac{k}{i}\right\rfloor * i k − ⌊ i k ⌋ ∗ i
进一步,原式= ∑ i = 1 n k − ⌊ k i ⌋ ∗ i = \sum\limits_{i=1}^n{k - \left\lfloor\frac{k}{i}\right\rfloor * i} = i = 1 ∑ n k − ⌊ i k ⌋ ∗ i 其中k k k 是个常数,可以拿出来,因此继续化简得到n ∗ k − ∑ i = 1 n ⌊ k i ⌋ ∗ i n*k - \sum\limits_{i=1}^n{\left\lfloor\frac{k}{i}\right\rfloor * i} n ∗ k − i = 1 ∑ n ⌊ i k ⌋ ∗ i .距离数列分块的形式还多了一个i i i ,对于求和的每一个单独的元素来说,i i i 就是一段等差数列求和的值.对于[ l , r ] [l,r] [ l , r ] 这段区间,i i i 部分产生的贡献即( l + r ) ∗ ( r − l + 1 ) / 2 (l +r) * (r - l + 1) / 2 ( l + r ) ∗ ( r − l + 1 ) / 2 再加上数论分块的部分即可.
代码
typedef long long ll;
#define int ll
int n,k;
int f(int n)
{
int res = 0;
for(int l = 1,r;l <= n;l = r + 1)
{
r = k / l ? min(k / (k / l),n) : n;
res += (k / l) * (r + l) * (r - l + 1) / 2;
}
return res;
}
signed main()
{
scanf("%lld%lld",&n,&k);
int fres = f(n);
printf("%lld",n * k - fres);
return 0;
}